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ABSTRACT 

We present measurements of the longitudinal magnetic field in the circumstellar environ- 
ment of seven classical T Tauri stars. The measurements are based on high-resolution circular 
spectropolarimetry of the He 1 A5876 emission line, which is thought to form in accretion 
streams controlled by a stellar magnetosphere. We detect magnetic fields in BP Tau, DF Tau 
and DN Tau, and detect statistically significant fields in GM Aur and RW Aur A at one epoch 
but not at others. We detect no field for DG Tau and GG Tau, with the caveat that these objects 
were observed at one epoch only. Our measurements for BP Tau and DF Tau are consistent, 
both in terms of sign and magnitude, with previous studies, suggesting that the characteris- 
tics of T Tauri magnetospheres are persistent over several years. We observed the magnetic 
field of BP Tau to decline monotonically over three nights, and have detected a peak field 
of 4kG in this object, the highest magnetic field yet observed in a T Tauri star. We combine 
our observations with results from the literature in order to perform a statistical analysis of 
the magnetospheric fields in BP Tau and DF Tau. Assuming a dipolar field, we determine a 
polar field of ^ 3 kG and a dipole offset of 40° for BP Tau, while DF Tau's field is consistent 
with a polar field of ~ — 4.5kG and a dipole offset of 10°. We conclude that many classical 
T Tauri stars have circumstellar magnetic fields that are both strong enough and sufficiently 
globally-ordered to sustain large-scale magnetospheric accretion flows. 

Key words: accretion, accretion discs - stars: circumstellar matter - stars; magnetic fields - 
stars: pre-main-sequence 



1 INTRODUCTION 

It has been postulated that accretion in low-mass pre-main- 
sequence stars is magnetica lly-controlled (e.g. Koenigl .199 1.; 
ICollier Cameron & CamDbeiii ri993). with the magnetic field dis- 
rupting a Keplerian disc and the accreting material plummeting 
along the field lines onto the stellar surface. This model provides 
an attractive solution to the angular momentum dissipation that is 
required to produce a slowly rotating protostar, and synthetic emis- 
sion line profiles produced assuming a magnetospheric accretion 
via a dipolar field are able to reproduce some characteristics of 
the o bserved profiles (e.g. Hartmann et al. 1994; Muz erolle et aL 
Il998h. 



One of the key requirements of this model is the presence 
of a strong, relatively structured, magnetic field a few stellar radii 
above the star. Traditional methods for measuring surface magnetic 
fields rely on the Zeeman effect. Actual splitting of absorption line 
profiles is rarely observed, because other broadening mechanisms 
(rotation, pressure, turbulence etc) usually dominate over the rel- 
atively weak Zeeman effect. Nonetheless, if the intrinsic shape of 
the profile is well understood, it is possible to attribute the resid- 
ual broadening to Zeeman splitting, and hence measure a mag- 
netic field strength. This method has been used with some suc- 
cess on classical T Tauri stars (CTTSs), and surface fields of a few 



kilogauss have been determined iBasri et al]|l992l : lGuenther et alJ 
[1999; Johns-Kiull et al. 199§). 

Although surface magnetic fields have been measured for a 
few CTTSs, it is not clear how the strength of these fields changes 
with radius and how ordered they are (e.g. Safiet 1998). Measure- 
ments of circumstellar magnetic fields require observations of the 
Zeeman effect in an emission line, and here the Zeeman broadening 
technique fails - primarily because the lines are broadened by the 
bulk motion of the accreting gas, and the intrinsic profile of the line 
cannot be determined with any certainty. Fortunately Zeeman split- 
ting also results in a polarization signature, with the cr-components 
of the transition being circularly polarized in opposite senses. The 
separation between these compo nents is directly proportional to the 
mean longitudinal field strength iBabcoclJl947il . 

Johnstone & Pension ( 1986, hereafter JP) attempted to mea- 
sure the magnetic fields in three T Tauri stars using this effect, 
and found a 2.3-a detection (— 548 G) for RU Lup, which wa s 
not confirmed by later observat i ons jjohnstone & Penstoiilll987h . 
More recently 'Johns- Krull et alj il999t) measured the Zeeman shift 
in circular polarization of the He I A5876 in BP Tau, using an 
echelle spectrograph on a 2.7-m telescope at McDonald Obser- 
vatory. The He I A5876 line is thought to form above hot impact 
regions on t he stellar surface in streams of gas being accreted by 
the star te.g. lValenti et alll993l) . lEdwards et aljjl994 found that 
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the line showed an inverse P Cygni profile for stars with high lev- 
els of continuum veiling - two well-known accretion indicators. 
IJohns-KruU et alJ jl999h inferred a mean longitudinal magnetic 
field of 2.4 kG, which was the first direct evidence of magnetic ac- 
cretion in a CTTS. Similar observation s o n three more CTTSs w ere 
presented bv lJohns-Krull et alj j200lh and lValenti et all <2003ft . 

Here we present measurements of the longitudinal magnetic 
field in seven CTTSs using the He I A5876 emission line. The aim 
of these observations was to expand the sample of CTTSs with 
known fields and to examine the long-term (many rotation cycles) 
stability of the magnetosphere. 



He I emission lines <Edwards et ai]|l994 . The TEK4 CCD detec- 
tor was used during the 2001 run, providing 0.40 A per pixel re- 
ciprocal dispersion; the newer MARCONI2 CCD was employed in 
2002 to obtain 0.22A per pixel. The resolution achieved was 0.7A 
(FWHM of the arc lines). Maximum sensitivity of the system was 
selected by inserting a circular-polarizing filter during setup, and 
finding the quarterwave-plate angle that gave the maximum con- 
trast between the intensity of the o- and e-beams. Observations of 
the target objects (typically 1000-s per frame) were bracketed by 
exposures with copper-neon and copper-argon comparison lamps 
to provide a wavelength calibration. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Zeeman effect 

In this paper we exploit the Zeeman effect to measure the mag- 
netic field around CTTSs. If there is a net magnetic field along the 
observer's line of sight, a spectral line is split and the separation 
in wavelength between the two cr-components is a measure of the 
field strength. The wavelength shift is 



AA 



■ 9.34 X lO^^A^ffoffB^mAkG"^ 



(1) 
(2) 



where A is the rest wavelength of line (in A), gcs is the effec- 
tive Lande g-factor of the line transition and is the net longi- 
tudinal magnetic field (Mathvs 1991). The value of g^s adopted 
for He I A5876 was 1.11, but the line is formed by multiple 
transitions and Qcs is dependent on the environment (see e.g. 
|johns-KrulletalJll999h . For other spectral lines, Lande g-factors 
were obtained from the solar line list table of I.F. Donati repro- 
duced at http : //bass2000 .bagn . obs-mip . f r. 



2.2 Observations 

Spectra were obtained with the ISIS spectrograph at the Cassegrain 
focus of the 4.2-m William Herschel Telesco pe. The configuration 
for cir cular spectropolarimetry is described in lTinbergen & RuttenI 
il992h . A quarterwave plate first converts from circular to linear 
polarization. Light from the target object then passes through one 
aperture of a comb dekker; the other apertures sample the back- 
ground sky emission. A comb dekker with 20-arcsec apertures was 
used for all exposures except the tungsten lamp flat-fields. The star 
and sky beams all pass through a Savart plate which spatially sepa- 
rates the orthogonal linear polarizations into o- and e-beams, each 
of which then passes through the standard spectrograph optics and 
is recorded on a CCD detector. Rotating the quarterwave plate be- 
tween exposures, by 90° relative to the first position, swaps the 
polarization sense of each beam and should reverse direction of the 
measured wavelength shift. Many of the possible systematic errors 
in the optics or detector would lead to unequal inferred magnetic 
field magnitudes from a pair of exposures, and this is the first of 
several tests we can use to assess the veracity of our results. 

The observations were made on nights beginning 2001 De- 
cember 3-5 and 2002 December 21-23 with the R1200 grating 
on the ISIS red arm. Approximately half of each run was lost to 
weather or technical problems. We observed calibration targets, 
and seven bright CTTSs in the Taurus- Auriga region with known 



2.3 Data reduction 

The procedure for obtaining a magnetic field measurement requires 
cross-correlation of two spectra recorded simultaneously on differ- 
ent areas of the CCD. Differential errors in the extraction or cali- 
bration of the spectra will create a false magnetic field detection; it 
is important to test the systematic effects are well below the (less 
than one pixel) expected Zeeman wavelength shift. 

Each CCD frame was first bias-subtracted and flat-fielded. The 
flat-field frames had a large intensity gradient towards the edges 
where vignetting was significant, so these frames were normalised 
so that they had a value close to 1 across the central area, where the 
spectral lines of interest were recorded. In the 2002 run, flat-field 
frames were recorded at both quarterwave-plate rotation positions 
(So, 6o +90°) and the appropriate master flat-field used for correct- 
ing each of the target exposures in a pair. Similarly, the arc com- 
parison frames were recorded at both quarterwave-plate angles, and 
the appropriate arc spectrum was used for calibrating each target 
spectra. 

The 2001 observations contained calibration frames made at 
only one rotation angle of the quarter-wave plate, and it was neces- 
sary to test the effect of reducing data from both quarter-wave plate 
positions with a flat-field from just one position. Test reductions 
were carried out with the 2002 data using mismatching flat-fields 
and the derived wavelength shifts were not significantly different 
from the results using flat-fields from both quarter-wave plate po- 
sitions. These results were also unaffected when a single arc-lamp 
exposure was used to calibrate each complementary pair of target 
exposures, so we had confidence in using these techniques for the 
2001 data. 

For each exposure, the two stellar spectra (corresponding to 
the left- and right-handed polarisation states) were extracted, as 
well as two spectra sampling the sky background. The sky spectra 
were of very much lower intensity than the stellar continuum and 
were unstructured over the wavelength extent of the He I A5876 
line. Tests on sample frames showed that sky subtraction did not 
affect the results; it was decided that this process only contributed 
noise and should be omitted from the data reduction sequence. 

Spectra of the arc lamps were obtained from their frames us- 
ing the same extraction windows used for the stellar spectra, then 
a third-order polynomial wavelength scale was fitted using iden- 
tifiable emission lines as references. The typical RMS error from 
the fit to the arc lines was ~ O.OSA. The four stellar spectra (an o- 
and an e-beam recorded at each quarterwave plate position) were 
rebinned to a common linear wavelength scale in preparation for 
determining the wavelength shift. The scale was chosen to be sim- 
ilar to that of the individual spectra. 
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2.4 Cross-correlation analysis 

A range of wavelength bins containing the emission or absorp- 
tion line of interest were copied from the full spectra for cross- 
correlation using the procedure described in JP. We manually in- 
spected the spectra to choose the limits for the line region, rather 
than employ an automated procedure. The emission lines in the T 
Tauri-class spectra were generally strong compared to the contin- 
uum, so most of the photons were recorded in a narrow wavelength 
range and it was not felt necessary to include the line wings where 
it was difficult to identify the transition to the local continuum in- 
tensity. 

The two rebinned spectra from each CCD frame were cross- 
correlated by evaluating a statistic, S, for one spectrum being 
used as a model for the other (after allowing for intensity rescal- 
ing). As described in JP, we made whole-pixel shifts between the 
spectra and fitted a parabola to the three values of S that were clos- 
est to the minimum value, S(min). The minimum turning point of 
the parabola gives the wavelength shift between the two spectra, 
AA, that is used to obtain the magnetic field strength (Equa- 
tion|2}. TableQlists for each of the target exposures. 

Although we had propagated uncertainty estimates for each 
pixel or wavelength bin throughout the data reduction, when 
the spectra were rebinned to a common wavelength scale it be- 
came non-trivial to estimate the uncertainty contribution from each 
source bin to each target bin because the errors are then corre- 
lated. We_fliereforeusedthe ratio of variances method (described 
in JP jLampton et aljl97dh to estimate the l-cr uncertainty interval 
for each Zeeman wavelength shift. The interval is found from the 
cross-correlation shift statistic, S, relative to S(min), its value at 
the minimum turning point of the parabola, such that: 



S < S(min) 



1 + 



N' 



-F(p, N' 



P, aj 



(3) 



where p is the number of free parameters used to find S(min); 
N' is the number of independent wavelength bins in the spectra; a 
is the confidence level required (i.e. ~ 0.68 for the l-cr interval); 
and -F(p, N' — p, a) is the Fischer-Snedecor ratio of variances 
function. One spectrum from each pair was scaled in intensity to 
match the other, and then shifted in wavelength, so p = 2. We 
follow the advice of JP, who show that after rebinning the spectra, 
the effective number of independent wavelength bins is reduced 
to 5/6 of the number of bins in the raw spectra. The interval in 
Equation |3| defines a range of wavelength shifts that gives (with 
Equation |2} the \-a uncertainties on the B^ values presented in 
TableQ 

A hybrid scheme was also tried which used inverse variance 
weighting when calculating the x? statistics. The advantage of this 
alternative was that the variance associated with each wavelength 
bin included information about the flat-field uncertainties, detector 
noise etc. The results were not significantly altered by this change 
in procedure, and we concluded that the original cross-correlation 
method was adequate. 

We found that small changes made to our procedures (such as 
altering the wavelength limits of the emission line, or the trials us- 
ing single calibration frames in Section l23l typically perturbed the 
individual measurements by approximately lOOG. When the un- 
certainties in the polarization standards (Section|23} are also taken 
into account, we estimate that the systematic errors in our results 
are probably of magnitude ~ 200G and we would not confidently 
claim to detect magnetic fields of lesser strength. We have not at- 
tempted to account for the systematic effects when quoting uncer- 



tainty intervals in this paper; they are always the formal errors de- 
rived from the cross-correlation procedure. 

The magnetic field of a target T Tauri star was expected to 
vary little over the duration of the observations (less than ~ Ih) 
on a given night. We therefore obtained a nightly Bz measurement 
for each star by computing a mean value weighted by the inverse 
variance of the individual estimates. FigQshows both the individ- 
ual measurements and the mean values, with their l-cr uncertainty 
intervals. We present the combined results in Table |2| along with 
characteristic stellar data from the literature. 

To test our estimates of the 1-sigma uncertanties we can calcu- 
late a reduced chi-squared statistic for all our measurements, with 
respect to the mean for the object on the night in question. This 
yields a xl of 0.836 with 59 de grees of freedom. Since we expect 
81 per cent of such experiments to exceed this value, it suggests 
that both our assumption that the polarization is unchanging within 
the ~lh of our observations, and our calculation of the uncertain- 
ties is broadly correct. Table|2|also gives the values of xi for each 
object on each night. Although there are a couple of uncomfortably 
high values of for individual objects on individual nights, given 
the above result, and the low number of degrees of freedom for the 
nightly chi-squareds, their statistical significance is unclear. Fur- 
thermore, the highest value corresponds to the data with the small- 
est uncertainty (86G) and thus we would expect it to be affected 
by the systematic uncertainties that we estimate to be about 200G. 
We therefore conclude that the x^ analysis supports our conclusion 
that there are systematic effects at the 200G level, but above this 
level our uncertainties are probably correct. 

2.5 Test stars 

We observed three targets with dependable levels of circular po- 
larization to test our instrumental setup and data reduction pro- 
cess. The Sun has a magnetic field strength much smaller than 
the sensitivity of our equipment, so it was used as a null refer- 
ence. Sunlight reflected from Vesta provided practical access to the 
solar spectrum; for the purposes of this experiment, the asteroid 
can be considered to be a passive reflector of circularly-polarized 
light. In poor seeing conditions, only two measurements had rea- 
sonable signal-to-noise ratios and these two useful exposures were 
obtained with the quarterwave plate at the same angle, so we can 
not test that the same magnitude of wavelength shift is measured 
as the beams are reversed (as described in Section l231 . An inverse- 
variance weighted mean of the shift in the Na I 5889A line repre- 
sented a mean field strength = — 20± 180G, which is consistent 
with the expected null result. 

Observations of the magnetic Ap star 53 Cam permitted com- 
parison with measurements published by other authors. Data ob- 
tained in bo th the 200 1 and 2002 runs were phased using the 
ephemeris of iHill et alJ il99 8). Our two measurements are shown 
in Fig. |2| along with those from two other sources. Rather than 
choose a sign convention for B^ based on the spectrograph setup, 
we have assumed that our 53 Cam measurements are consistent 
with the previously published data, and should therefore both have 
a positive sign. Our data then are of the magnitude expected from 
the previous studies, indicating that our measurement technique is 
reliable. The sign convention adopted is later shown to also be con- 
sistent with previous measurements of T Tauri magnetic fields pub- 
lished by other authors. 

Observations were made of a linear polarization standard star 
to determine whether our instrumentation would apparently detect 
a Zeeman wavelength shift because of unwanted conversion from 



4 N.H. Symington, T.J. Harries, R. Kurosawa and T. Naylor 



Figure 1. Plots of the (G) results from individual exposures. The circles are the measurements derived from individual exposures; the horizontal lines are 
the 1-(T uncertainty intervals for those values. Each night's data for each object were combined using an inverse-variance weighted mean. The vertical solid 
and dashed lines show the mean and l-cr intervals respectively. 
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linear and circular polarization in the optics. The star HD2 5443 
(5.1 per cent V band linear polarization: Tumshek et alll99(J) was 
observed during the 2001 run and several spectral lines were cross- 
correlated. The wavelength shift found by our data reduction pro- 
cess would be equivalent to a magnetic field of — 194±45G, which 
would be a significant detection. It is possible that this star does 
have an intrinsic longitudinal magnetic field, but it is more likely 
that the apparent wavelength shift between the two beams is caused 
by imperfect polarization optics. We expect that any effects of po- 
larization cross-talk are going to be within the 200G systematic 



error estimate we discussed above and do not threaten the validity 
of our experimental results. 



3 RESULTS 

Here we describe our measurements of individual objects, and re- 
view the observational evidence for magnetically controlled phe- 
nomena. 
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Table 1. Summary of T Tauri star obsei-vations and individual results. Times quoted are UT at exposure midpoint. Quarterwave plate rotation angles (0) are 
quoted in degrees, relative to the adopted zero point (Section [2.2r Uncertainties in the longitudinal magnetic field, Bz, are the \-a interval. 
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Figure 2. Magnetic field measurements of the Ap star 53 Ca m. Our data are 
sho wn as filled squares , along with those o et aljfl993) (open squares) 
and lWade et'ai]l2OO0h (open circles). 



3.2 Possible detections 

In 2001 we detected a magnetic field on RW Aur (+0.46±0.21 kG) 
although it was at the limit of our setup's sensitivity. Subsequent 
observations on two nights in 2002 found B^, close to zero. On 
the night of 2002 December 22, for GM Aur was found to be 
— 1.26 ± 0.46 kG - a significant detection - but by the following 
night was consistent with zero (+0.49 ± 0.42 kG). 

These stars mig ht be expected t o sho w evidence of mag- 
netospheric accretion. iGullbring et alj il998h found the accretion 
rate of GM Aur to be of the order of ~ 1O~®M0 yr~^and 
iMuzerolleetail fcOOl') presented fits to the Balmer lines and NaD 
profiles of RW Aur based in a magnetospheric accretion model with 
M =~ 10~* Mq yr~^. High resolution time-series spectroscopy 
of RW Aur A revealed periodic modulation of the emission line 
profiles (P = 2.77 A, Petrov et al. 2001). Non-axisymmetric accre- 
tion, due to either a close companion or an offset dipole configura- 
tion, was postulated to explain the varia bility. Magnet ic fields are 
also implicated in the formation of its jet iHirth et alll9 94). 

We are less confident in our claims for these CTTSs than for 
others (Section f3.3> . The behaviour may be evidence for rotational 
modulation of the magnetic fields, but detections at more than one 
epoch would be desirable for these two targets. 



3.1 Null detections 

Two of the target CTTSs had results that were consistent with zero 
longitudinal field. Our cross-correlation analysis (Section IZ?! leads 
to formal 2-a upper limits of IB^j < 0.80 kG for DG Tau and 
\Bz\ < 0.32 kG for GG Tau. Each of these stars was observed 
on only one night, so it is possible that the measurements were 
made at epochs when the orientations of the magnetic fields were 
unfavourable. 



3.3 CTTS with detected magnetic fields 

|johns-KrulletalJ<1999h used circular spectropolarimetry to search 
for magnetic fields on BP Tau. They determined an upper-limit 
of 200 G for the surface field, but found = 2460 ± 120 G in 
the magnetosphere from the He I A5876 line. They concluded that 
the accretion must be confined to streams with small footp rints on 
the stellar surface. Further He! A5876 observations I Valenti et alj 
2003) showed that B^ is variable on a time-scale of days, indicating 
rotational modulation of a structured magnetosphere. 



6 N.H. Symington, T.J. Harries, R. Kurosawa and T. Naylor 



Table 2. Nightly re sults for each star. The v sini data are from Hartmann & Stauffer '1989). Other properties for each object were from references: 
l.|a)hen & Kuhl ll979ft : 2.fGunbrine et^ ll998ft : 3. Hartigan et al. 1 1995); 4. Valenti et al. 1 1993). Mass accretion rate (M) is in units of IQ-^Mq yr-\ 
V sin j in kms"'^. 
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1.57 










2002 Dec 21 


6 


3820 


480 


0.106 










2002 Dec 22 


4 


1980 


290 


0.361 










2002 Dec 23 


6 


1170 


150 


0.093 


DFTau 


M0.5(l) 


16 


1.767 (2) 


2002 Dec 21 


6 


-980 


160 


0.121 










2002 Dec 23 


6 


-717 


86 


2.40 


DG Tau 


K7-M0 (3) 


22 


20 (3) 


2002 Dec 21 


6 


420 


570 


0.151 


DN Tau 


MO 


10.2 


0.033(2) 


2002 Dec 23 


6 


630 


160 


2.20 


GG Tau 


K7-M0(l) 


10.2 


0.175(2) 


2002 Dec 23 


5 


70 


250 


0.795 


GM Aur 


K7-M0(l) 


12.4 


0.096(2) 


2002 Dec 22 


5 


-1260 


460 


0.825 










2002 Dec 23 


4 


490 


420 


1.98 


RW Aur A 


K3 (4) 


19.8 


3.4 (4) 


2001 Dec 05 


4 


460 


210 


0.263 










2002 Dec 21 


3 


210 


960 


0.271 










2002 Dec 23 


6 


190 


330 


0.258 



Our BP Tau observati ons from 2001 and 2002 confirm those 
of I.Tohns-Krull et alJ h999l) and I Valenti etail OoOS), although our 
peak Bz, at 4kG, is the strongest yet found for any CTTS. In our 
sign convention (see Section |231 . Bz was always positive, as it 
was in the previously published results. The field strength appeared 
to be monotonically declining over the three consecutive nights of 
the 2002 observing run. Three data points are not suffi cient to draw 
any c onclusions about periodicities in the data, but the lValenti et alJ 
i2003 ') observations also showed a smooth change over six consec- 
utive nights (with a minimum value of 0.3 kG), so it likely that the 
dominant magnetic time-scale is many days, possibly matching the 
stellar rotation period (7.6 d; Simon et al. 1990). We investigate this 
possibility further in Section l?!4l 

DF Tau is a binary system ( Chen et al although the pri- 

mary i s responsible for the continuum and line emission in the blue 
iLamzin e t al. 2001). I Valenti et alj 12003") observed a longitudinal 
field strength that varied smoothly between — . 3 and —1.0 kG over 
six nights. Our two measurements are within this range, so we have 
another indication that the T Tauri magnetospheres have persistent 
long-term characteristics. 

It should be made clear that Fig.Qappears to show the indi- 
vidual DF Tau measurements from 2002 December 23 alternating 
between two separate magnetic field strengths. The bimodal dis- 
tribution is accounted for by the change of quarterwave-plate an- 
gle between exposures. The sign of the magnetic field, but not its 
magnitude are expected to change after that procedure. Small-scale 
magnitude variations were seen in other targets, but the later set of 
DF Tau observations shows the most significant effect. It is likely 
that imperfections in the spectrograph optics are affecting either 
the polarization signature of the star, or introducing a systematic 
error into our wavelength calibrations. The mean magnitude of Bz 
in these data is sufficiently high that we are confident in ruling out 
a null detection, but the formal error estimates in the final result 
should be treated sceptically. 

Our measurements of Bz for DN Tau show an overall statis- 
tically significant detection (+0.63 ± 0.16kG) although several of 
the contributing individual data points were consistent with zero 
magnetic field. We suggest that we have probably made the first 
detection of a longitudinal magnetic field from this star. 

Two of the CTTSs with detected mag netic fields have 
shown evidence for surface hotspots: DN Tau JVrba et aljfigsS 



Table 3. Published longitudinal magnetic field measurements for the 
T Tauri stars BP Tau an d DF Tau using the He I A5876 line. References: 



l. lValentietall<2003h : 2. |johns-Krull et alJ<199Sft . 



Name 


Night 


Bz(G) 


± l-a 


BPTau 


1997 Nov 21 


2460 


120 




1998 Nov 26 


2740 


170 




1998 Nov 27 


1490 


170 




1998 Nov 28 


850 


150 




1998 Nov 29 


330 


280 




1998 Nov 30 


390 


320 




1998 Nov 31 


610 


150 


DFTau 


1998 Nov 26 


-300 


110 




1998 Nov 27 


-520 


90 




1998 Nov 28 


-810 


80 




1998 Nov 29 


-940 


90 




1998 Nov 30 


-1020 


330 




1998 Nov 31 


-1000 


110 



Ref. 



iBouvier et all 198^ and DF Tau jBouvier et alJ 1993ft - which has 
been mapped bv lunruh et all il998l) using Doppler tomographv. 



3.4 Models for BP Tau and DF Tau 

Our observations of the longitudinal magnetic field of T Tauri stars, 
have some targets in co mm on with those recently published by 
I.Iohns-KrulI et alJ h 999ft and I Valenti eTall i2003ft . Those authors' 
results for the stars BP Tau and DF Tau are shown in Tabl^ With 
all three data sets, there is a sufficient number of measurements (11 
and 8 respectively of BP Tau and DF Tau) to permit a preliminary 
statistical analysis of the magnetic field properties. 

We investigated the case where the measured field arises 
solely from an accretion stream following the field lines of a dipo- 
lar magnetosphere. The line-of-sight component of such a magnetic 
configuration, Bz, is given by: 



Bz = B 



pole 



'cos i COS P + sin i sin f3 cos t 



(4) 



where Bpoic is the field strength at the poles of the magneto- 
sphere; i is the observer's inclination with respect to the rotation 
axis; f3 is the dipole offset - the inclination of the magnetic polar 
axis with respect to the rotation axis; and (j) is the rotational phase 
angle at the time of observation. 



T Tauri stellar magnetic fields: He I measurements 
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Figure 3. Kolmogorov-Smimov statistics for BP Tau (left) and DF Tau (right). Each pixel shows the probabihty that the observations are consistent with a 
simulated distribution of measurements from a dipole magnetic field with given values of -Bpoic ^nd dipole offset (/3). Black represents a probability of > 50 
per cent, white is 0, and the greyscale shows a linear range between these two values. The step changes seen in the BP Tau simulation are a consequence of 
employing the K-S test with a small number of data points. 



A Monte Carlo computer code was used to make simulated 
observations at many phases based on input parameters of Bpoic, 
i and /3. The relative spacing (but not actual values) of the phase 
points were chosen to reflect the timing of the real observations, 
i.e. a random phase point was first chosen to represent 'night 1' of 
an observing run, and then the elapsed time to the next observation, 
divided by the rotation period, was used t o derive the subseq uent 
phase points. The periods used were 7.6 d (Simon et al. 199(1), for 
BP Tau and 8.5 d for DF Tau (Bouvier 1990). The inclination pa- 
rameter, i, was allowed to vary by 10° about fixed values from the 
literature: 50° for BP Tau, 70° for DF Tau iValenti et all2003h . 

Many such sets of observations were simulated for varying 
values of Bpoic and /3 and a Kolmogorov-Smimov (K-S) test was 
carried out to determine whether the genuine observations were 
consistent with the simulated distribution. In Fig|3| the K-S proba- 
bility is shown as i?poie and /3 are changed, and it can be seen that 
the models are ruled out for many combinations of the these param- 
eters. In each case, there is a practical constraint that Bpoic must be 
at least as great in magnitude as the maximum observed value of 
■ Imposing a narrow range of viewing inclinations will also in- 
troduce a degeneracy between Bpoic and j3, i.e. a Bpoic value that 
does not match the star's actual -Bpoic might be compensated for 
by changing the dipole offset so that the distribution of simulated 
observations is more widely spread and then is consistent with the 
real observations. 

The BP Tau observations are found to be best matched by 
Bpoic ~ 3kG and P ~ 40°. DF Tau was consistent with ~ 
— 4.5kG and (3 <~ 10° - the difference in dipole offset is con- 
sistent with our expectation, since our model gives smaller dipole 
offsets for objects with smaller variability for obvious physical rea- 



4 CONCLUSIONS 

Measurements were made of the mean longitudinal magnetic field, 
Bz, from Zeeman splitting of the He I A5876 line in the spec- 
tra of seven classical T Tauri stars. For the star that had been 
previously best-studied with this technique, BP Tau, the field 
strengths w ere mostly within the range found by previous authors 
jjohns-KruU et al.. 1999:.Valenti et aL.2003i) . but we also measured 
the strongest field so far detected for any T Tauri star (4kG). An- 
other star, DF Tau, was found to have a magnetic field that was also 
within the range of previous measurements. A first Bz detection 
was made for DN Tau, and the stars GM Aur and RW Aur each 
showed one measurement that was above our instrumental sensitiv- 
ity limit. In total, five of the seven targets showed evidence for a 
magnetic field at one or more epochs. 

There are now a total of nine CTTSs in the literature with 
precise magnetospheric field measurements, and of this sample a 
statistically significant field has been detected in seven objects. In- 
terestingly our repeat observations of BP Tau and DF Tau indicate 
that the fields may have characteristics that are persistent on a time- 
scale of years (i.e. hundreds of rotation cycles). It is evident that 
classical T Tauri stars have fields with the magnitude {r^ 1 kG) and 
stability expected in theoretical studies of the magnetospheric ac- 
cretion process (e.g. Romanova et al. 2003). 

We note that the dipole offset estimates derived for BP Tau 
and DF Tau (40° and < 10°respectively) would result in very dif- 
ferent magnetospheric strea m structures acco r ding to the magneto- 
hydrodynamical models of iRomanova et all i2003h . The simula- 
tions suggest that DF Tau should accrete from two diametrically 
opposed streams, while BP Tau (with its larger dipole offset) 
would have a more complex accretion structure. However, it ap- 
pears BP Tau sho ws line profile variabili ty that is consistent with 
an inclined dipole jGullbring et'S]ll99d) . while the variability of 
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DF Tau is apparently diffic ult to interpret in terms of th e magneto- 
spheric accretion paradigm Jjohns-KruU & Basrll99'i . 

The next stage of these investigations should be a simultane- 
ous measurement of the mean field and the longitudinal field, per- 
haps using time-series measurements to map the surface magnetic 
field structure using a technique such as Zeeman Doppler Imaging 
(ZDI: Semel 1989.) . 
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